Ha Emission Variability in the 7-ray Binary LS I +61 303 
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ABSTRACT 

LS I +61 303 is an exceptionally rare example of a high mass X-ray binary (HMXB) that also 
exhibits MeV-TeV emission, making it one of only a handful of "7-ray binaries" . Here we present 
Ha spectra that show strong variability during the 26.5 day orbital period and over decadal time 
scales. We detect evidence of a spiral density wave in the Be circumstellar disk over part of 
the orbit. The Ha line profile also exhibits a dramatic emission burst shortly before apastron, 
observed as a redshifted shoulder in the line profile, as the compact source moves almost directly 
away from the observer. We investigate several possible origins for this red shoulder, including an 
accretion disk, mass transfer stream, and a compact pulsar wind nebula that forms via a shock 
between the Be star's wind and the relativistic pulsar wind. 



Subject headings: accretion, accretion disks - stars: winds, outflows - stars: emission line, Be 
individual(LS I +61 303) 



stars: 



1. Introduction 

LS I +61 303 is a high mass X-ray binary 
(HMXB) that is also a confirmed source of very 
high energy 7-ray emission. The system con- 
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sists of an optical star with spectral type BO Ve 
and an unknown compact c ompanion in a highly 
eccentric, 26.5 day orbi t (lAragona et all 12009: 



Grundstrom et al.l 120071: ICasares et alTl2005h 



While the system has a relatively low X-ray 
luminosity for a HMXB, LS I +61 303 is the 
15th brightest 7-ray source included i n the Fermi 



LAT 1-year Point Source Catalogue (jAbdo et al 
20101) . The Be disk interacts with the compact 
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companion, producing orbital phase modulated 
emiss ion across the electrom agneti c spectrum: 



TeV (lAlbert et all 120061 12008ft. GeV (lAbdo et al. 
20091). X-ra v dParedes et al.l Il997t fLeahvl l200lh " 
optical Ha (Grundstrom et al.l 120071). and radi o 



(jGregorv fc Tavlorlll978l;lTavlor fc Gregory||1982D . 
Periodic radio emission from LS I +61 303 was 



first reported bv lTavlor fc Gregory! (|1982l ). They 
defined the arbitrary reference for zero phase at 
HJD 2,443,366.775 that remains the conventional 
definition for this system. The binary orbital 
period (P = 26.4960 ± 0.0028 d) has been de- 
termined from the periodic radi o outbursts tha t 
peak near </>(TG) = 0.6 - 0.8 (|Gregorvl Eool . 
and this period is confirmed by optical spec- 



troscopy (jAragona et al.l 120091) . Periast ron occurs 
at 6(TG) = 275 (lAragona et al.ll2009h . 



A number of authors have observed an over- 
all increase in the Ha emission line strength and 



the phase of radio maximum (Paredes et al.l 1994 


Zamanov et al. 1999; Liu et al. 2000t Liu & Yan 


2005 


). 


Grundstrom et al. 


2007: IStovanov et al. 


2008 


Paredes et al. ( 1994 


) and Liu et al.l (2000) 



proposed that an accretion disk around the 
compact companion may contribute to the Ha 
emission at this phase. On the other hand, 
Grundstrom et al.l (|2007l ) proposed that this pe- 
riodic emission feature is attributed to a spiral 
density wave in the Be star disk that is accreted 
by the compact companion shortly after perias- 
tron. There is a correlation between the X-ray 
flux and the disk emissi on strength during the or - 



bit (see Figure 5 from iGrundstrom et al.l [2007) 



The large orbit-to-orbit variations in Ha are quite 
similar to the variability of the radio outbursts. 

We recently obtained an extensive collection of 
red optical spectra of LS I +61 303 to determine 
an upd ated orbital ephemeri s for the spectroscopic 
binary (jAragona et al. 2009). These observations, 
reviewed in §2, also recorded the evolution of the 
Ha emission during a full orbital cycle. In §3 we 
discuss the short-term (intra-orbit) and long-term 
variability of the Ha emission. Finally, we discuss 
the origin of the emission in §4. 

2. Observations 

During 2008 October and November, we ob- 
tained 83 red optical spectra of LS I +61 303 at 



the KPNO coude feed (CF) telescope over 35 con- 
secutive nigh ts. These observation s are described 
in detail by lAragona et al ] (|2009h . To summa- 
rize, the spectra covered a wavelength range of 
6400-7050 A with a resolving power R = A/AA ~ 
12, 000. We generally obtained 2-3 spectra of LS 
I +61 303 each night. The spectra included only 
two prominent stellar features, Ha and He I A6678, 
as well as the interstellar line at 6613 A. The CF 
spectra were calibrated using standard procedures 
in IRAP0 and were interpolated onto a common 
heliocentric wavelength grid. To rectify the spec- 
tra to a unit continuum, we omitted a wide region 
surrounding the emission line features to avoid any 
possible contamination of the emission line profile 
shapes. 

We observed LS I +61 303 again on UT dates 
2009 October 17-19, or 0.654 < 0(TG) < 0.733, 
at the Wyoming Infrared Observatory (WIRO) 
using the Long Slit Spectrograph. We used the 
LS-1 grating in 2 nd order, with the KG3 order- 
sorting filter, and these spectra cover the wave- 
length range 5400-6800 A with R ~ 4500. Due 
to variable clouds each night, we used exposure 
times of 300-900 s to obtain a total of 10 spectra 
during the run. We also obtained a CuAr lamp 
exposure just before each science exposure. The 
WIRO spectra were bias corrected, flat fielded, 
and wavelength calibrated using standard proce- 
dures in IRAF. Finally, we rectified the WIRO 
spectra to a unit continuum using line-free regions 
and interpolated them to a common heliocentric 
wavelength grid. 

3. Ha Emission Properties 

Since th e orbital period of LS I +61 303 is only 
26.4960 d (|Gregorvll2002l) . our 2008 collection of 
CF spectra enable the first detailed study of the 
Ha emission throughout an entire, single orbit. 
We show in Figure [1] the Ha line profiles and a 
gray-scale image of this line over the span of that 
observing run. Some glitches that do not affect the 
emission line profiles have been removed for clar- 
ity. Note that neither the line profiles nor gray- 
scale plots are folded by orbital phase, but rather 
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they reveal true chronological variations in the line 
profile behavior as a function of H JD and the cor- 
responding orbital phase. 

A prominent emission feature stands out as a 
"red shoulder" in the Ha line profile between pe- 
riastron and apastron. The feature first appears 
in our spectra from UT date 2008 November 6, at 
0(TG) = 0.592. By the following night, the red 
component is broader still. The extension sub- 
sides over the next two nights, and vanishes by 
UT 2008 November 10, 0(TG) = 0.740. Although 
this unusual red shoulder was detected in eight 
spectra over four nights, we were initially con- 
cerned that the spectra might be contaminated 
by scattered light or some other instrumental ef- 
fect, or that it could be an artifact of our rectifi- 
cation procedure. We examined all of the raw CF 
spectra carefully, including other targets observed 
on the same nights using the same instrumental 
configuration, and we repeated our rectification 
of the complete set of LS I +61 303 spectra to 
exercise maximum caution in the continuum fit- 
ting process. In fact, this unusual red shoulder 
is not the result of contamination or a transient 
optical flare; it is a persistent fea t ure near this or- 
bital phase dParedes et al.lll994t |Liu et al.l l2000t 
iGrundstrom et alj|2007h . 

The Ha line profiles of the WIRO spectra are 
qualitatively similar to the 2008 CF spectra at 
comparable orbital phases, shown in Figures [5] 
and [3] A similar, but significantly weaker, ex- 
tended red-shifted profile is marginally visible near 
<HTG) ~ 0.6-0.7 during 1998 August, 1999 Octo- 
ber (immediately followed by the nearly complete 
ionization of the disk), 1999 Novem ber, and 2000 
December ( Grundstrom et al 1 l2007h . It is appar- 
ent that over the past decade, the orbit-to-orbit 
variations in the Ha emission line profile arc quite 
large. 

We measured the equivalent width of Ha, Wn a , 
for each spectrum by directly integrating over the 
line profile. (We use the convention that Who is 
negative for an emission line.) The errors in Wn a 
are typically about 10% due to noise and place- 
ment of the continuum. Figure |4] shows that dur- 
ing our Coude Feed run, Wn a decreases slightly 
just before periastron. Since Wna is correlated 
to the radius of a Be st ar's circumstellar disk 
( Grundstrom fc Gies 20061 ). we interpret the de- 
cline in emission as a slight decrease in disk radius 



as gas is stripped away by the compact compan- 
ion. Wn a then rises dramatically with the on- 
set of the red shoulder emission component near 
4>(TG) ~ 0.6. 

Figure [4] also com pares our recent W n a wit h 
those measured by IGrundstrom et al.l ([20071 ) . 



Their data were accumulated over six different 
observing runs over 1998-2000, and the long term 
differences in emission strength are substantial. 
Most of their runs do not provide an opportunity 
to compare the emission strength at both the pe- 
riastron and subsequent quadrature phases. How- 
ever, their first run over HJD 2451053-2451065 
(1998 August) has the best coverage of this or- 
bital range, and the stronger emission strength 
at periastron suggests that the Be star's disk was 
slightly larger in radius than in 2008. The sub- 
sequent rise in the 1998 emission occurs sooner 
than in 2008, but gaps in the 1998 coverage do 
not allow us to determine whether the emission 
peak was comparable. 

Stovanov et al.l ( 20081 ) also obtained 110 Wn a 



measurements over multiple years (see their Fig- 
ure 2). Although they do not distinguish between 
orbital cycles, they find that that the emission can 
peak anywhere between 0.35 < <£(TG) < 0.7. The 
strong emission episodes never seem to occur past 
apastron. 

Smaller temporal changes in the 2008 CF spec- 
tra suggest additional Ha emission variability, 
so we subtracted the mean emission line profile 
to investigate the residu als carefully. Although 



Grundstrom et al.1 (|2007l ) did not observe any ra- 



dial velocity variations in the peaks or wings of 
Ha, we do see them in our dataset, especially 
in the central trough of the line. Therefore we 
formed a mean Ha line profile by shifting each 
spectrum to the star's rest velocity. Then we 
shifted the mean to the orbital velocity of the 
star a ccording to the ephemeris of lAragona et al ' 



(2009) and subtracted it to produce Ha emission 
residuals, or difference spectra, shown in Figure 
[5j We also produced a version of Figure [5] with- 
out accounting for orbital motion, and the results 
were not significantly different. During about half 
of the orbit, 0.9 < </>(TG) < 0.6, the difference 
spectra reveal a partial S-shaped pattern similar 
to a spiral density wa ve that is commonly ob - 
served in Be star disk s ( Porter fc Riviniusl 2003 ). 
Zam anov et al ] (Il99fil ) also observed a strong blue 
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peak near </>(TG) = 0.23, which supports the de- 
velopment of a spiral density wave near periastron. 
After this phase, the peculiar red shoulder devel- 
ops. 

The observed behavior of the spiral density 
wave and the red shoulder emission components 
are quite different. To investigate whether these 
features form in the circumstellar disk of the opti- 
cal star or in an accretion disk around the compact 
object, we used Gaussian fits of the peak residual 
emission in our difference spectra to measure its 
radial velocity, V r , and full width half maximum, 
FWHM. In some cases, the difference spectra did 
not clearly show a distinctive peak, so we omitted 
those spectra from measurement. We also omitted 
measurements in cases where the residual emission 
profile deviated substantially from a symmetric, 
Gaussian profile. Many of the red shoulder emis- 
sion profiles were not included as a result. 

Between the orbital phases 0.59 < 4>(TG) < 
0.63, our Gaussian fits of the red shoulder mea- 
sure its V r between 300-443 km s _1 . The op- 
tical star has V r ~ —52 near these phases, and 
a center of mass ve locity 7 = —41.41 km s _1 
( Aragona et al. 20091). Assu ming a 12.5M Q opti- 
cal star ( Casares et al.|[2005h and a 1.4M neutron 
star, the expected radial velocity of the neutron 
star (and any associated accretion disk) should 
only be ~ 54 km s _1 . If the V r of the red shoul- 
der was due to orbital motion, this would require 
a mass ratio q = M1/M2 > 25, which is unlikely. 
We conclude that the observed V r of the red emis- 
sion shoulder does not correspond to the orbital 
motion of the co mpact companion as proposed by 
Liu et all (120001 ). 

The FWHM of the residual emission peaks are 
shown in Figure [6l We estimate that the errors in 
our FWHM measurements are about 50 km s _1 
due to noise intrinsic to the difference spectra and 
minor deviat ions from the Gau s sian p rofile used 
in our fits. iGrundstrom et all ( 20071 ) measured 
a projected rotational velocity Vs'mi = 104 ± 5 
km s _1 for LS I +61 303. The true equatorial 
rotational velocity is V eq , and the base of the 
Keplerian circumstellar disk should rotate at the 
critical velocity V C rit, where Vea/Vcrit ~ 0.85 for 



phases < 0(TG) < 0.5 and 0.7 < 0(TG) < 1.0, 
the low FWHM is generally consistent with a spi- 
ral density wave in the Be star disk. On the 
other hand, the FWHM of the red shoulder emis- 
sion component is significantly higher, suggesting 
a more turbulent region of gas. 

The Who and FWHM of the red shoulder emis- 
sion component suggest that this orbital feature 
does not originate from within the Be star's disk. 
From the V r of this feature, we also rule out orbital 
motion of an accretion disk surrounding the com- 
pact companion. In the next section, we discuss 
other possible origins for this feature. 

4. Discussion 



Paredes et al.1 (|1994l) first proposed that red 
emission may be due to an ac- 
around the compact companion. 



shoulder of 
cretion disk 

Romero et all ([20071 ) present a three-dimensional 
smoothed particle hydrodynamics simulation of 
accretion onto the neutron star in LS I +61 303. 
They predict that spiral structure is induced both 
in the accretion disk and the Be circumstellar 
disk near periastron. In their simulations, the ac- 
cretion rate peaks with a phase delay of about 
0.3P or b- The short duration of the red shoulder 
emission occurs about 0.3 phase past periastron, 
consistent with their simulations. However, Ha 
emission from an accretion disk is not generally 
visible in the optical spectra of other HMXBs, and 
the radial velocity of the red shoulder is inconsis- 
tent with the orbital motion of an accretion disk, 
so the accretion disk itself is an unlikely origin for 
this emission feature. 



a typical Be star (jMcSwain et all 12008V There- 



fore emission that originates in the circumstel- 
lar disk of LS I +61 303 should have FWHM 
< 2V sin i/(Veq/V cr u) = 329. During orbital 



Ra ther than an accretion disk source. IGrundstrom et al 
(|2007l ) propose the development of a tidal stream 
within the Be star disk near periastron, which ex- 
tends beyond the truncation radius of the Be star 
disk and into the vicinity of the compact compan- 
ion. Such a t idal stream is consis tent with the 
simulations of Romero et all ( 2007 ). When the 
red shoulder appears at 0(TG) of 0.59-0.73, we 
may be looking somewhat down the throat of the 
resulting mass stream as the compact companion 
is receding from the observer; the relative orbital 
geometry is shown i n Figure [7] using the recent 
orbital ephemeris by I Aragona et al I (l2Q09h . 

If the red shoulder of emission is produced by a 
tidal stream, we should observe a correlation be- 
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tween the disk radius just before periastron (as 
measured from Whq at this phase) and the emis- 
sion strength of the mass stream. This is not ob- 
served over multiple epochs of Wh q measurements 
(Fig. 0]). Furthermore, the stream should be rel- 
atively collimated as it flows across the inner La- 
grangian point. The dynamical trajectory of the 
stream should not vary from orbit to orbit, so the 
timing of the emission peak should be constant 
during each orbit, but there is no evidence of any 
consistency over many orbits. While a tube of 
gas would appear optically thick while it is ori- 
ented along our line of sight, the emitting area of 
the mass stream is unlikely to be large enough to 
cause up to 50% increase in the emission strength 
relative to the normal Be star's disk emission. Fur- 
thermore, the broad FWHM of this feature does 
not support a collimated mass flow onto the neu- 
tron star (Fig. [6|). 

The red shoulder emission could be formed as 
the induced tidal stream falls back onto the Be 
circumstellar disk. Since the more massive Be 
star dominates the gravitational field, most of the 
tidally disturbed material will not be accreted, in- 
stead crashing back onto the Be disk. Around 
<fi(TG) ~ 0.6, the infalling gas would be moving 
away from the observer and would impact the red- 
shifted part of the Be disk. The density would 
increase quickly as the plume hit the disk and 
the emission would increase because of the den- 
sity squared dependence. If this origin is correct, 
the orbit-to-orbit variation of the red shoulder be- 
havior may be due to the changes in the Be disk 
gas density at the azimuth and time near the po- 
sition of the companion at periastron. This may 
be quite variable since the gas involved would be 
in the tenuous, outer reaches of the disk. If the 
red shoulder does come from the outer part of the 
Be disk facing the companion, the mean densities 
there may be too low to produce the He I 6678 
emission (usually formed in the inner, high den- 
sity regions of the disk). This is probably consis- 
tent with the lack of a red shoulder in He I A6678 
(jAragona et al.ll2009t Fig. 7). 

We can estimate the infall velocity of the tidal 
stream by assuming that its circular velocity is 
comparable to the periastron velocity of the neu- 
tron star, 281 km s _1 . Disk gas would be launched 
outwards at periastron, and then fall back towards 
the other side of the disk with a radial velocity of 



about +281 — 41 = +240 km s 1 (correcting for 



the s ystemic velocity of the binary jAragona et al 



2009). Our first measurement of the red shoulder 
is at V r = +293 km s _1 (Table 1), which is larger, 
but the disagreement is not too bad considering 
our simple model. 

A fourth option for the origin of the Ha shoul- 
der could be a cone-shaped wind shock region. LS 
I +61 303 may contain a shrouded pulsar whose 
relativistic wind interacts with the optical star's 
wind to form a wind shock region (|Dubusl [20061: 



Dubus et al. 2010l ). Isolated pulsar wind neb- 



ulae occasionally exhibit Ha spectra typical of 
a Balmer-dominated shock (BDS), including the 
presence of both narrow (~ 10 km s _1 ) and broad 
(~ 500— 1000 km s" 1 ) hydrogen emission line com- 
ponents (Heng 2010h . While the blending of the 
Be star's emission with the red shoulder makes 
it difficult to identify both Ha components of the 
BDS in LS I +61 303, the unusual broadness of the 
red shoulder emission is consistent with a BDS. 
While many pulsar wind nebulae do not exhibit 
such a BDS, the high density of the stellar wind 
surrounding the pulsar wind nebula in LS I +61 
303 increase s the chances of forming such a shock 
(|Hendl2010h . The temporary nature of the red 
shoulder may suggest that the BDS only forms 
when the high density tidal stream interacts with 
the neutron star. Based on the large, positive V r 
of the red shoulder emission observed in this sys- 
tem, the flow of gas along the shock cap should 
be oriented away fr om the observer. However, 
Rom ero et al. I (|2007h modeled such a wind shock 
region in LS I +61 303, but they found difficulties 
reconciling the observed morphology of the radio 
emission and the spectral energy distribution with 
their models. 

Dubus et"aTl(|2010h advocate a cometary shock 



flow in LS I +61 303 trailing the neutron star's 
orbit (and thus opening toward the observer near 
0(TG) ~ 0.6) if the compact object interacts with 
gas from the Be circumstellar disk that is mov- 
ing more slowly in a corotating frame. According 
to their interpretation, we should observe blue- 
shifted emission from the gas flow, contradicting 
the observed red shoulder. The orbital velocity of 
the neutron star is given by 



'orbit — 



GA1* 



1 



(M* + M NS ) 2 \tns (ins, 



(1) 
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where Mms and are the masses of the neutron 
star and optical star, tns is the neutron star's dis- 
tance from the center of m ass, and clns is the neu- 
tron star semimajor axis ([Hilditchl 120011 ). Com- 
pare this to the approximately circular Keplerian 
orbit of the circumstcllar disk material at the neu- 
tron star's location, 



isk 



GM* 



These velocities are equivalent for 



r NS = 2a NS 



M+ + M NS 



o-ns 



(2) 



(3) 



Assuming M* = 12.5M Q (jCasares et al.ll2005l) and 
M NS = 1.4 M Q with P = 26.4960 d, then a N s = 
80.9 Rq. The neutron star would move faster than 
the circumstellar disk for 0.194 < (f>{TG) < 0.356, 
only very close to periastron. Accounting for only 
the circular component of the neutron star's or- 
bit (rather than the total orbital speed) changes 
these results only slightly. If the Keplerian disk 
extended out to the neutron star's location, the 
disk would actually stream past the slow moving 
neutron star near apastron, potentially reversing 
the orientation of the shock flow in agreement with 
our observations of the Ha emission. However, the 
Be circumstellar disk should be truncated within 
the periastron distance of the neutron star. Any 
stellar wind interacting with the neutron star near 
(j)(TG) ~ 0.6 would be oriented radially away from 
the Be star, so any comet-shaped interaction re- 
gion would be skewed away from the line of sight. 

If the red shoulder forms in an optically thin 
region, we should observe flux from other orbital 
phases unless the emission is only temporary. A 
short-duration red shoulder might be formed by 
any one of several mechanisms: the formation of a 
short-duration emitting structu re, like a bullet o f 
gas being ejected (as in SS 433; iGies et al.l 120021) . 
a strong beaming effect that causes the feature 
to be observed only at particular orientations, or 
the complete ionization of the emitting structure 
at most other orbital phases. Since we do not 
observe relativistic velocities, a bullet of ejected 
gas or beaming effects are unlikely. Ionization ef- 
fects would make more sense if the emission was 
observed more symmetrically about apastron. In 
our opinion, the most likely mechanism to form a 



short-duration emission structure occurs when the 
induced spiral dens ity wave extends acr oss the full 
binary separation (jRomero et al.1 120071 ). Whether 
the tidal stream falls back onto the Be circumstel- 
lar disk or interacts to form a temporary pulsar 
wind nebula is somewhat inconclusive. We recom- 
mend further simulations to study the evolution of 
the tidal stream throughout the entire orbit, and 
we urge further observations in an effo rt to con- 



clusiv ely detect a pulsar in this system (|Rea et al 
l2010h . 
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Table 1 
Properties of Ha Emission 



HJD Orbital Wrq FWHM of Difference Spectrum V r of Difference Spectrum 



(-2,450,000) 


Phase (TG) 


(A) 


(km s" 1 ) 


(km s" 1 ) 


4756.7256 


0.87 


-9.14 


86.0 


223.3 


4756.7482 


0.87 


-9.90 






4757.6644 


0.91 


-11.06 






4757.6855 


0.91 


-9.52 






4758.6607 


0.94 


-9.89 


90.3 


179.9 


4758.6827 


0.94 


-10.22 


87.1 


183.6 


4759.6877 


0.98 


-10.48 






4759.7088 


0.98 


-10.20 






4760.6803 


0.02 


-10.64 






4760.7014 


0.02 


-10.62 






4761.7508 


0.06 


-10.88 


166.7 


78.0 


4761.7719 


0.06 


-10.74 


151.6 


74.9 


4762.6955 


0.10 


-11.59 






4762.7166 


0.10 


-11.45 






4763.7018 


0.13 


-11.49 


177.6 


-180.5 


4763.7229 


0.13 


-11.12 


169.0 


-182.3 


4764.6373 


0.17 


-11.18 






4764.6583 


0.17 


-11.46 


157.2 


-219.3 


4765.0282 


0.18 


-11.01 


192.7 


-246.7 


4765.6837 


0.21 


-11.02 


187.3 


-281.9 


4765.9010 


0.22 


-11.22 


144.3 


-304.3 


4765.9221 


0.22 


-10.97 


149.7 


-308.9 


4766.7893 


0.25 


-10.54 


133.5 


-339.9 


4766.8574 


0.25 


-10.31 


140.0 


-339.0 


4766.9291 


0.26 


-10.60 


143.2 


-339.5 


4767.7007 


0.28 


-10.21 






4767.7805 


0.29 


-10.60 






4767.8459 


0.29 


-10.91 


80.8 


-398.4 


4768.7456 


0.32 


-10.56 






4768.8149 


0.33 


-10.78 






4768.8841 


0.33 


-10.57 






4769.7368 


0.36 


-10.59 


139.9 


-69.5 


4769.7970 


0.36 


-10.39 


102.2 


-59.5 


4769.8738 


0.37 


-9.82 


136.6 


-65.8 


4770.7309 


0.40 


-11.19 


158.1 


-34.3 


4770.7911 


0.40 


-11.24 


160.3 


-36.6 


4770.8624 


0.40 


-11.29 


150.6 


-32.1 


4771.8542 


0.44 


-11.13 


150.6 


-32.1 


4771.8753 


0.44 


-11.00 


111.9 


-45.7 


4773.7536 


0.51 


-12.02 


104.3 


190.9 
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Table 1 — Continued 



HJD 

(-2,450,000) 


Orbital 
Phase (TG) 


(A) 


FWHM ot Difference Spectrum 
(km s ) 


V r of Difference bpcctrum 
(km s ) 


4773 7752 


0.51 


— 12.41 


107.5 


193.6 


4774 Q024 


56 

yj .OK) 


— 12.71 


1 33 3 


181 7 


4775 7528 


59 


— 1 6 47 


362 1 


293.2 


4775 775Q 


0.59 


— 16.02 


329.9 


295.5 


477fi 7648 


0.63 


— 16.35 


471.6 


424.8 


477fi 7862 


0.63 


— 16.08 


498.4 


435.7 


4777.7764 


0.67 


— 14.55 






4777 8Q1 Q 


0.67 


— 14.76 






4778 7466 


0.70 


— 13.08 


108.6 


— 15.1 


4778 7677 


0.70 


— 13.08 


100.0 


—8.7 


4779 6636 


0.74 


— 10.42 






4779 6851 


0.74 


— 10.56 






4781 71 32 


81 


—9 39 






4781 7344 

'-k I O J. • 1 04:4: 


81 

U.Ol 


—9 74 






4782 6531 


0.85 


—9.28 


116.2 


—86.4 


4782 6747 


0.85 


—9.30 


135.6 


—80.5 


4782 8258 


0.86 


—8.74 






4782 8470 


86 

u.ou 


—8 84 

O.Ot: 






4784 6Q52 


93 


— 1 37 


79 6 


— 17 


4784 71 63 


0.93 


—9.91 






4784 7536 


0.93 


— 10.09 






4784 7747 


93 

yj .iJO 


— 10.11 






4785 8022 


97 


—9 03 






4785 8233 


97 

yj . iJ i 


—9 1 7 






4786 7401 


00 


—9 70 






4786 761 2 


0.00 


—9.63 








n rizL 




on a 


1 f\f\ ft 

— ±00. o 


4787.7129 


0.04 


-9.82 


88.3 


-163.2 


4787.7340 


0.04 


-9.56 






4787.7551 


0.04 


-9.32 


107.6 


-169.1 


4787.7762 


0.04 


-9.71 


91.5 


-175.5 


4788.8069 


0.08 


-10.94 


106.6 


-200.1 


4788.8280 


0.08 


-11.00 


86.1 


-204.3 


4789.6443 


0.11 


-10.93 


139.9 


-215.7 


4789.6654 


0.11 


-10.88 


120.6 


-212.5 


4790.7113 


0.15 


-9.90 


146.4 


-238.1 


4790.7324 


0.15 


-10.26 


145.3 


-239.0 


4790.7536 


0.15 


-9.95 


162.6 


-227.1 


4790.7747 


0.16 


-9.94 


165.8 


-238.1 


4791.7059 


0.19 


-10.58 


180.9 


-276.9 
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Table 1 — Continued 



HJD 


Orbital 


w Ba 


FWHM of Difference Spectrum 


V r of Difference Spectrum 


(-2,450,000) 


Phase (TG) 


(A) 


(km s- 1 ) 


(km s- 1 ) 


4791.7270 


0.19 


-10.27 


148.6 


-286.5 


5121.8113 


0.65 


-14.80 






5121.8591 


0.66 


-14.27 






5121.8812 


0.66 


-14.35 






5122.8172 


0.69 


-14.60 






5122.8282 


0.69 


-14.64 






5122.8562 


0.69 


-14.43 






5122.8644 


0.69 


-14.60 






5122.8900 


0.69 


-14.46 






5122.8978 


0.70 


-14.45 






5123.9006 


0.73 


-14.51 
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4750 7 




VELOCITY (km s 1 ) 



Fig. 1. — The upper plot shows the Ha line pro- 
file of LS I +61 303 over our continuous 35 nights 
of observation, sorted by HJD, and the lower plot 
shows a gray-scale image of the same line. Note 
that the lower plot of the gray-scale spectra are 
not folded by orbital phase but are placed in the 
same chronological order as the upper plot, with 
the orbital phases indicated. The intensity at each 
velocity in the gray-scale image is assigned one of 
16 gray levels based on its value between the min- 
imum (dark) and maximum (bright) observed val- 
ues. The intensity between observed spectra is cal- 
culated by a linear interpolation between the clos- 
est observed phases. The solid wh ite line shows 



the t heoretical V r curve solution (jAragona et al 
20091) . 




Fig. 2. — The mean CF Ha line profiles from 
<j)(TG) = 0.59 (dotted line), 0(TG) = 0.63 (dot- 
dashed line), <f>(TG) = 0.67 (dashed line), and 
4>(TG) = 0.70 (dot-dot-dot-dashed line) are com- 
pared to the mean CF spectrum (solid line). All 
spectra have been shifted to the rest velocity of 
the optical star. 
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- 1 f— 




0.5 



6520 



_i i_ 



6540 6560 6580 

WAVELENGTH (A) 



6600 



Fig. 3. — The mean WIRO Ha line profiles from 
</>(TG) = 0.66 (dotted line), 4>{TG) = 0.69 (dot- 
dashed line), and </>(TG) = 0.73 (dashed line) arc 
compared to the mean CF spectrum (solid line). 
All spectra have been shifted to the rest velocity 
of the optical star. 
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Fig. 4. — Variations in Wu a function of or- 
bital phase. Measurements from the six runs of 
Grundstrom et al. (2007; identified by the date 
HJD-2,450,000) and the mean Wn a from each 
night of our 2009 WIRO run are shown folded 
by orbital phase. The consecutive set of 2008 CF 
measurements have not been folded but are placed 
in their true chronological order. The times of 
periastron and apastron are marked with vertical 
dotted lines. 
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VELOCITY (km s 1 ) 



Fig. 5. — Emission residuals, or difference spectra, 
in the same format as Figure Q] The difference 
spectra have been smoothed for clarity. 
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Fig. 6. — FWHM of the peak residual emission 
measured from Gaussian fits of the difference spec- 
tra. 
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Fig. 7.- 
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O rbital geometry of L S I +61 303, 
adapted from Aragona et al. (2009). The relative 
orbits (r/a) of the 12. 5 Mm optical star (dotted 



line; Casares et al. 2005h and a presumed 1.4M Q 
neutron star companion (solid line) are shown. 
The center of mass is indicated with a cross, and 
the thin solid line is the orbital major axis. Points 
around the orbit indicate steps of 0.1 in orbital 
phase, and the phases of periastron and apastron 
are labeled. The stars move counterclockwise in 
the figure. 
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